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Abstract 


Ab initio calculations at the self-consistent-field and singles plus doubles 
configuration-interaction level are used to determine accurate spectroscopic param- 
eters (D e , r e , u e ) for most of the alkali and alkaline-earth fluorides, chlorides, 
oxides, sulfides, hydroxides and isocyanides. Numerical Hartree-Fock (NHF) calcu- 
lations are performed on selected systems to ensure that the extended Slater basis 
sets employed for the diatomic systems are near the Hartree-Fock limit. Extended 
gaussian basis sets of at least triple-zeta plus double polarization quality are 
employed for the triatomic systems. By dissociating to the ionic limits, most of the 
differential correlation effects can be embedded in the accurate experimental elec- 
tron affinities and ionization potentials. With this model, correlation effects are 
relatively small (0.0-0.3 eV), but invariably increase D 0 . The importance of corre- 
lating the electrons on both the anion and the metal is discussed. 

The theoretical dissociation energies (D 0 ) are critically compared with the 
literature to rule out disparate experimental values. The theoretical studies com- 
bined with the experimental literature allow us to recommend D c values that are 
accurate to 0.1 eV for all systems considered. The systematic treatment of many 
different systems reveal many trends. For example, the dissociation energies of the 
alkali and alkaline-earth hydroxides are observed to be less than the corresponding 
fluorides by just slightly less than the difference in electron affinities of F and OH. 
In general, there is a strong correlation between the dissociation energy (to ions) 
and r e , because the bonding is predominantly electrostatic in origin. 

Theoretical 2 II— 2 E + energy separations are presented for the alkali oxides and 
sulfides. The ground states of all the alkali sulfides are shown to be X 2 II. An 
extensive study of the 2 II— 2 S + energy separation in KO reveals a 2 E + ground state 
at all levels of theory. The separation is shown to be sensitive to basis set quality, 
and in the NHF limit the 2 E + state is lower by about 250 cm -1 . The separation is 
almost unaffected when the 16 valence electrons are correlated at the singles plus 
doubles level using an extended Slater basis. 


PRECEDING PAGE BLANK NOT FILMED 
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I. Introduction 


The dissociation energies (D 0 ) of many of the ionic molecules containing the 
alkali and alkaline-earth atoms have been determined experimentally. The dissoci- 
ation energies of the alkali halides reported by Brewer and Brackett [l] in 1961 are 
still accepted today. Their dissociation energies are limited in accuracy only by the 
quality of the initial vapor pressure data and by the necessary corrections for gase- 
ous imperfection. In contrast, experimental dissociation energies for the alkaline- 
earth oxides show a wide variation. For MgO, the experimental D c values [2-7] 
range from 3. 5-4. 3 eV, whereas a recent theoretical study [8] obtains 2.65±0.16 eV. 
For CaO, two chemiluminescence studies yielded disparate D 0 values of 
(^4.76±0.15 eV [9] and 4.11±0.07 eV [10]. Theoretical studies [ll] gave excellent 
agreement with the lower value, supporting Dagdigian’s conclusion that the higher 
value was incorrect as a result of interference with CaCl [10], 


It is demonstrated in this paper that theoretical dissociation energies can be 
computed with 0.1 eV accuracy for large classes of ionic diatomic and triatomic 
systems. Considered herein are most of the alkali and alkaline-earth fluorides, 
chlorides, oxides, sulfides, hydroxides and isocyanides. The theoretical dissociation 
energies are capable of ruling out disparate experimental values, and allow an 
overall assessment of the various experimental methods used to determine dissocia- 
tion energies. Our theoretical methods described in detail in the next Section 
reflect the fact that the charge distribution in ionic molecular systems much more 
closely resembles the constituent ions than the neutrals. Hence, by dissociating to 
the ionic limits and then correcting to the neutral limits using the accurate experi- 
mental ionization potentials and electron affinities, the problem can be formulated 
such that relativistic and correlation contributions to the dissociation energies are 
relatively small (^0.3 eV). This is documented for a large number of systems by 
comparing self-consistent-field (SCF) and singles plus doubles configuration- 
interaction, CI(SD), values for D 0 . 


For the alkali oxides and sulfides we have considered both the lowest 2 S + and 
2 n states. This is of particular interest for the alkali oxides, where there is a 
change in ground-state symmetry from LiO( 2 II) to CsO( 2 £ + ). For the KO 
molecule there is both conflicting experimental [12-13] and theoretical [14-15] 
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evidence for the ground-state symmetry. This work provides rather strong support 
for KO(X 2 £~). In contrast, the alkali sulfides Li-Rb are shown to have X 2 II 
ground states, probably as a result of the increased bond lengths. 
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I. Methods 


The fraction (f) of ionic character in a bond is difficult to quantify. Pauling’s 
criterion [16], f=/x/er e , where n is the dipole moment and r e is the equilibrium 
internuclear separation is not a good measure of ionicity in highly ionic systems 
owing to the large deformations caused by the electrostatic field of the ions. 
Although Mulliken populations can show ionic character, they are difficult to 
quantify and often vary significantly with the one-particle basis (e. g. diffuse func- 
tions increase the overlap population which is difficult to assign). The best meas- 
ure of ionic character is probably the ratio r x /r e where r x is the hypothetical cross- 
ing point where a purely ionic potential curve crosses the asymptote of the 
covalent curve. In Table I we have tabulated this ratio for selected alkali and 
alkaline-earth fluorides, chlorides, oxides and sulfides. Systems for which r x /r e ^ 
2.5 potentially have some covalent character in the bond. By this criterion, all of 
the alkali and alkaline-earth oxides and sulfides as well as the alkaline-earth 
halides of Be and Mg could contain some covalent character. 

For highly ionic systems, such as the alkali halides, the molecules can be con- 
sidered to be composed of ions, each of which is polarized by the electrostatic field 
of the other. For these systems most of the binding energy arises from purely elec- 
trostatic interactions between the ions: charge-charge interaction, charge-dipole 
interaction, dipole-dipole interaction, and quasi-elastic energy stored in the induced 
dipoles. According to the Rittner model [17], this portion of the binding (<f>) can be 
expressed in terms of the atomic polarizabilities (a^ and as 


4 > 



e 2 (<*i+a 2 ) 

2r 4 


2e 2 a 1 Q 2 


(1) 


The Rittner model has yielded reliable binding energies for the alkali halides where 
accurate experimental r e values are available. In addition, the model of Kim and 
Gordon [18] for describing forces between closed-shell atoms has been successfully 
applied to the alkali halides. 

Probably the most rigorous theoretical treatment of the dissociation energies 
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(D e ) of selected alkali halide molecules is given in a series of papers by Matcha 

[19]. He computed D e at the SCF level with large Slater bases, and then corrected 

for correlation and relativistic effects using atomic data. In his formalism the total 

HF corr 

dissociation energy of the molecular system MX is given by D e = D e + D e + 
. D e re ^ where 


D corr p corr , p coit p corr p corr 
e ~ +&X 


( 2 ) 


and 


D/ el ~ E^ 1 -f-Ex el -E^- 1 (3) 

Formulated in this way Dj~ OTT is positive and can be quite large. Matcha 
obtained rather good agreement with the available thermochemical D 0 . His values 
are slightly low because his basis sets were still not at the HF limit, and he 
neglected the inter-fragment correlation energy (see later discussion). 

The formalism that we employ herein is essentially equivalent to Matcha’s, but 
avoids explicitly knowing the total correlation and relativistic energies of the 
atoms and corresponding ions. Also, our approach includes the interfragment 
correlation which, although small in total energy (e. g. 0.18 eV for NaF), is 
entirely differential since it vanishes in the asymptotic limit. The dissociation 
energy of MX can be written as 

D e = E(M + ) + E(X’) - E(MX,r e ) - IP(M) + EA(X) (4) 

Hence, we explicitly dissociate to the ionic limits and then correct to the neutral 
ground state atomic limits using the accurate experimental ionization potentials 
(IP) available from Moore [20], and electron affinities (EA) available from Hotop 
and Lineberger [21]. Since most of the differential correlation and relativistic 
effects are absorbed into the IP and EA, the SCF dissociation energies, D e , are 
potentially quite accurate. Note that it is necessary to apply this formalism to a 
molecular state that is well represented by the HF configuration. This is true of 
the ground states of all of the syst ems considered herein except for the alkaline- 
earth oxides and sulfides whose l £ + ground states require a multireference descrip- 
tion owing to the fact that they are a mixture of singly and doubly ionized struc- 
tures. Hence, for these systems we have applied our formalism to either the excited 
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a 3 n or A 1 IT states, both of which are well described by a single reference configura- 
tion. The ground-state D 0 can then be computed if the T g is known for the excited 
state. 


For the formalism in Eq. 4 to be accurate, the basis sets must approach HF 
quality. Since the atoms approach the HF limit much more quickly than the molec- 
ular systems, the D e values almost invariably increase monotonically with 
improvements in basis set quality. For all of the diatomics considered herein, we 
have employed extended Slater bases. They begin with the accurate sets of 
Clementi and Roetti [22] and were modified slightly where necessary to reduce 
problems with linear dependency. The basis sets are then further augmented with 
polarization functions to describe the considerable distortions that arise in the field 
of the ions. These basis sets are described briefly in Table II. Additional details can 
be found in forthcoming publications [23-28]. To assess how close our SCF energies 
are to the HF limit, we have considered two approaches. First, for the lighter 
molecular systems (e.g. NaF, NaO, KO and BeCl) we have carried out numerical 
Hartree-Fock (NHF) calculations using an implementation of McCullough’s numer- 
ical diatomic code [29] on the Cray XMP. The numerical calculations show that 
our Slater bases are very near the HF limit (error < 0.02 eV). For the heavier sys- 
tems such as SrO, we have used basis set saturation studies to assess the degree of 

basis set incompleteness. Here the error is somewhat larger (0.20 to 0.25 eV), but 

SCF . 


most of the loss is in the core, so that our D„ 


is only low by about 0.07 eV. 


Also, for the triatomic systems considered herein, we have used extended gaussian 

SCF 

basis sets, so again there is a tendency for our D e values to be slightly low. 


Since near HF quality D e values can be produced using the formalism in Eq. 4, 

HF 

the question naturally arises as to how well these D e reproduce the true values. 

HF 

In other words, the D e will be accurate if the two conditions 

E corr ^ p corr . p corr / c \ 

MX ~ ^M* +^X- l 5 ) 

and 

E rel ^ p r*l la\ 

MX ~ ^M'+^X- ( 6 j 

are satisfied. The second condition is expected to be well satisfied for all systems 
considered herein, since only the valence electrons are substantially distorted by 
the electrostatic field. We have examined the approximation in (5) by performing 
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singles plus double configuration-interaction, CI(SD), calculations from the SCF 
reference configuration for all systems as well as more restricted Cl calculations for 
selected systems. We find, invariably, that correlation increases the dissociation 
energies, making the D e lower bounds. Also, since E^®/ 1 is always quite small 
(0.0-0. 3 eV) the SCF model as used in Eq. (4) is a good approximation to D e . 
Apparently, SCF calculations using extended basis sets accurately portray the 
charge distribution in these ionic molecules near their equilibrium geometry. 

In all Cl calculations we correlate the n and (n-1) shells on the alkali and 
alkaline-earth cations, which includes the Is electrons for Li and Be. The outer- 
most seven valence electrons are correlated for the oxide and sulfide anions, and 
the outermost eight valence electrons are correlated for the fluoride, chloride, 
hydroxide and isocyanide anions. Hence, for example, 17 electrons are correlated 
for the alkaline-earth halides, except for Be where 11 electrons are correlated. 
When this many valence electrons are correlated, one must potentially contend 
with size consistency, basis set incompleteness, and basis set superposition errors, 
especially considering the small differential correlation contribution to D„. 

To facilitate the discussion of electron correlation let us divide the 
configuration-state functions (CSFs) into three classes. Class 1 and 2 contain all 
single and double excitations from the orbitals that can be identified as X" and 
M + , respectively. Class 3 contains the important pair-pair terms, which are the 
double excitations where one electron is excited from an orbital and the other 
from an X~ orbital. This third class is size consistent, since it contributes nothing 
at infinite separation, but increases in importance as the bond distance decreases. 

If only class 1 CSFs are included, the calculation is size consistent, but the 
energy is not invariant to a unitary transformation among the and X~ occu- 
pied orbitals. For example, accidental degeneracies between the M + and X* orbi- 
tals result in arbitrary mixings that tend to reduce the correlation energy of both 
the alkali oxides and fluorides. Hence when the electrons on only one center are 
correlated, we use corresponding orbitals [30] to rotate the orbitals to give max- 
imum overlap with the M~*" orbitals. If the orbitals are not rotated, then the effect 
of the class 1 CSFs, which is to decrease D 0 and u> e and to increase r e , can be signi- 
ficantly accentuated. If only the class 1 CSFs are included, then the resulting spec- 
troscopic constants (D 0 ,r e ,w e ) are in worse agreement with experiment than are 
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the HF values. These comments are illustrated by the Cl results for the 2 S + 
ground state of NaF in Table III. If only the eight valence electrons on F" are 
correlated. D e is reduced by 0.08 eV, and by an even larger amount if the orbitals 

o 

are not first, localized. Also, r e increases by 0.02A worsening the agreement with 
experiment. In contrast, correlating the eight valence electrons on Na has little 
effect on D e or r e , which is not surprising considering the compact nature of the 
Na"*” orbitals. 

The class 3 excitations (pair-pair terms) have an opposite and somewhat larger 
differential effect than the class 1 excitations. Hence, including both classes (still 
size consistent) tends to increase D e and w e slightly, and to decrease r e , resulting in 
better agreement with experiment, particularly for the heavier alkali and alkaline- 
earths. The sum of the three disjoint classes, which contains all of the CSFs in the 
full CI(SD), should approximate the full calculation. In practice, the sum is a 
slight over-estimation of the combined effect, because the importance of the pair- 
pair terms (class 3) is reduced somewhat in the full calculation. The fact that the 
effect of the three classes is nearly additive suggest that the calculations are nearly 
size consistent if compared with the supermolecule. One can also rationalize the 
absence of any significant size-consistency error based on the small differential 
correlation energy observed. This is further supported by the fact that corrections 
for higher excitiations using either the Davidson’s formula [31] or CEPA1 [32] have 
negligible effect on the computed D e . The energy of the supermolecule was based 
on a point computed at an intemuclear separation of at least 20 Bohr. Note that 
there is a rather large size-consistency error («0.5 eV) if the calculations are refer- 
enced to the sum of the ions. 

Another potential problem that arises in the Cl calculations is basis set super- 
position errors. At this level of correlation treatment, significant demands are put 
on the quality of the basis set. For example, it is necessary to have sufficiently 
tight polarization functions to correlate the contracted spatial distribution of the 
occupied orbitals on M~**- We have computed the basis set superposition errors in 
our large Cl calculations using the counterpoise method. This approach is an 
upper bound to the error, since the counterpoise method tends to overestimate the 
superposition error, especially if the occupied space of the ghost atom is not 
excluded. The absence of a single basis function (such as a compact 3d function on 
M) can result in several tenths of an eV superposition error, and artificially 
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increase the calculated D e . These statements are well illustrated by the following 
example for CaF. If we use three Slater 3d functions for Ca atom (a=4. 76666, 
2.1222, 0.8007) we Find essentially no basis set superposition error at the SCF 
level, but a 0.129 eV error at the CI(SD) level when the fluorine ghost basis is 
brought up to 3.7 Bohr. When the Slater 3d basis is expanded to four Slater func- 
tions (q=6.0, 3.0, 1.5, 0.75) this error is reduced to 0.029 eV, and further to 0.014 
eV when the occupied fluorine orbitals are excluded. Reducing the superposition 
error by 0.10 eV decreased the dissociation energy by 0.08 eV and increased r e by 

c 

0.012A . Using the counterpoise method, we observe a 0.035 eV superposition error 
for fluorine in the presence of the Ca basis (0.030 eV with the calcium occupied 
space deleted). Hence, with the larger CaF basis, the total superposition error is 
about 0.05 eV at the CI(SD) level. Based upon tests for several systems, we 
believe our basis set superposition errors are less than 0.1 eV and, to a large 
extent, are cancelled by basis set incompleteness. However, this example illustrates 
the necessity of using very complete basis sets at the CI(SD) level to keep basis set 
superposition errors below the small differential correlation effect on D e . 

In performing our CI(SD) calculations from the SCF reference configuration for 
open shells, we generally invoked the interacting space [33], i.e. only those double 
excitations that have non-zero matrix elements with the SCF reference configura- 
tion are included. This has a small effect on the total energy, but essentially no 
effect on the spectroscopic parameters. Although, in general, we included all vir- 
tual orbitals in the Cl, our results indicate that eliminating the high-energy virtu- 
als (core correlating orbitals) also has very little effect on the spectroscopic param- 
eters. We also find that the spectroscopic parameters are insensitive to corrections 
for higher excitations [31], which is not surprising, considering the small total Cl 
effect observed. We have performed CEPAl calculations [32] on several systems, 
but have not reported these results here since they are not significantly different 
from the CI(SD) values. 

The initial phase of the study was done on the Cray-XMP using the DERIC- 
SWEDEN or MOLECULE-SWEDEN [34] molecular structure codes. When the 
Cyber-205 became available at NASA-Ames, we began calculations using the 
Karlsruhe adaptation [35] of the COLUMBUS codes [36-39]. The Slater integral^ 
were again evaluated with the diatomics integral program DERIC [40]. Extensive 
tests were performed to ensure that the two independent sets of programs gave 
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identical CI(SD) energies. More details of the program comparisons are given in 
the original papers [23-28]. 
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III. Results and Discussion 


Before proceeding with a discussion of the dissociation energies of the indivi- 
dual molecular systems, it is useful to discuss the overall trends that have been 
observed. In Fig. 1 we have plotted the dissociation energies, D e , of the alkali and 
alkaline-earth fluorides, chlorides, oxides, sulfides and hydroxides with respect to 
ground state ions versus equilibrium bond length. The nearly linear relationship 
between bond length and D e is striking. There is a tendency of the molecular sys- 
tems containing Li and Be to fall above the curve, which probably arrises from 

either the tendency of these systems to have a larger component of covalent char- 

2 2 6 

acter or the different core structure on the metal, i.e. s versus s p . The similarity 
of the fluorides and hydroxides is also quite evident from Fig. 1. Although the 
correlation diagram is not sufficiently accurate to predict D e values, it is capable of 
ruling out widely divergent ones, such as all of the experimental values for MgO. 

In the following Sections where we discuss the different classes of molecular sys- 
tems, we present D 0 values at both the SCF and CI(SD) level. In almost every 
case correlation acts to slightly increase D 0 and decrease r e , generally improving 
agreement with experiment. As discussed previously, the D 0 are probably 
lower bounds, and the CI(SD) D 0 are more likely low than high, especially for the 
heavier metals where we do not fully account for the bond contraction. We feel 
that all of the CI(SD) D 0 values reported herein for the diatomic systems (except 
those involving Cs and Ba) are accurate to 0.1 eV. Our results on triatomic sys- 
tems (hydroxides, isocyanides) that employ extended gaussian basis sets may be 
slightly less accurate, but again are more likely too low than high. Our theoretical 
results in conjunction with the available experimental values allow us to recom- 
mend reliable D 0 values for all systems considered in this study. 

This study also reports values for the spectroscopic parameters r e and u e . The 
CI(SD) results are generally in excellent agreement with experiment, although 
there is a tendency for theory to obtain slightly long bond lengths for the heavier • 
systems. Note that there is several thousands of an Angstrom uncertainty in the 
theoretical r e values, and a much larger (~30 cm -1 ) uncertainty in the w e values 
owing to the procedure of Fitting the theoretical data to a parabola. In some 
instances we have done additional points on a tighter (0.05 Bohr) grid to stabilize 



- 13 - 


the values. 


A. Alkali fluorides 


Two complete sets of experimental D 0 values exist for the alkali fluorides. 

298 

Brewer and Brackett [l] deter min ed D 0 values using a third-law treatment of 
vapor pressure data combined with the enthalpies of formation of the solid alkali 
halides and the gaseous monatomic elements. Their experimental numbers are lim- 
ited in accuracy by the original vapor-pressure data, and by the necessary correc- 
tions for gaseous imperfection. For an accurate comparison with our theoretical 
values given in Table IV, one should subtract 0.02-0.03 eV from their values to 
correct to absolute zero. The second set of experimental D 0 values given in Table 
IV are the flame photometric values of Bulewicz et al [41]. These values are in sub- 
stantial agreement with the thermochemical values, but contain a large uncer- 
tainty due to the insensitivity of the flame photometric method for the fluorides. 
For NaF a chemiluminescent determination [42] led to a lower limit (D 0 °^ 5.33 
eV) in good agreement with the flame photometric result [41], but considerably 
larger than the generally accepted thermochemical value [l]. The bond energy of 
CsF is well established experimentally since the bond energies determined by ther- 
mochemical [l], photoionization [43], flame photometry [41] and collision-induced 
dissociation [44] techniques give essential agreement. 

Comparing our SCF and CI(SD) D 0 values in Table IV, we Find that correla- 
tion uniformly increases D„, and that this effect is somewhat larger for the heavier 
alkali atoms. We have not yet performed the CI(SD) calculation for CsF, but 
predict a value of about 5.23 eV. This is in excellent agreement with experiment, 
except for the older thermochemical value of Scheer and Fine [45], which would 
appear to be about 0.2 eV too low. Overall, our CI(SD) D 0 values are in excellent 
agreement with experiment, especially the thermochemical values of Brewer and 
Brackett [l]. 

The CI(SD) D c values in Table IV should be accurate to about 0.1 eV. Hence, 
we recommend adopting 5.0 ±0.1 eV for NaF, which supports the thermochemical 
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value and is within the error bars of the flame experiments, but which rules out 
the chemiluminescent value of Ham [42]. We conclude, therefore, that their exists a 
small (0.3 eV) energy barrier in the formation of the emitting state in his reaction. 

The theoretical r e and u; e values are compared with experiment in Table V. The 
agreement with the experimental r e and w e values is excellent. There is a tendency 
at the SCF level to overestimate the bond length for the heavier alkali atoms 
where the pair-pair correlation is more important. The fact that the CI(SD) calcu- 
lations do not fully correct the deficiencies of the SCF may be due to basis set limi- 
tations or to the CI(SD) procedure itself, which may underestimate slightly the 
important differential effect of the pair-pair terms. As one can see from the strong 
correlation of D e and r e in Fig. 1, it is likely that the theoretical D 0 values are 
slightly small for the heavier systems where the CI(SD) r e values are slightly larger 
than experiment. 


B. Alkaline-earth fluorides 


Our theoretical dissociation energies for the alkaline-earth fluorides (Be-Sr) are 
in excellent agreement with the available experimental values. For BeF we favor 
the older mass spectrometric determination of Hildenbrand and Murad [46] over 
the later determination by Farber and Srivastava [47]. Our D 0 for BeF includes 
correlating the Be Is electrons, which produces a D 0 almost 0.1 eV greater than 
calculations in which the Is is not correlated. For MgF our D 0 of 4.66 eV lies 
between the two mass spectrometric determinations [48-49]. Since further basis set 
saturation is likely to increase our value slightly, we favor a value of 4.7±0.1 eV, 
which is consistent with all available values. 

Our CI(SD) D 0 for CaF lies within the error bars of the old thermochemical 
value of Blue et al. [50-51], but clearly suggests that the flame photometry D 0 of 
Ryabova and Gurvich [52] is too large. Apparently, the D 0 values for the fluorides 
determined by flame photometry are systematically large (except for LiF), but are 
usually correct within the rather large error bounds given. We would recommend 
adopting our theoretical value or a value slightly larger for CaF. For SrF, our 
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CI(SD) D 0 value agrees with all three experiments within their stated error bounds 
if one allows a 0.1 eV uncertainty in our value. However, our results again suggest 
that the flame photometry value [52] is too large and that the Ehlert et al. [48] 
mass spectrometric determination is too low. We are in excellent agreement with 
the mass spectrometric determination of Hildenbrand (D o =5.62±0.07) [49], and 
would recommend adopting this value for SrF. Work is presently in progress 
[23,53] on BaF. Preliminary results are in excellent agreement with experiment 
[48,49,51]. 


C. Alkali chlorides 


For the alkali chlorides rather reliable experimental values are available for 
comparison. The flame photometry technique [41] is much more sensitive for the 
chlorides, and produces values in good agreement with the thermochemical values 
[1]. Although our CI(SD) D 0 values in Table VI are in excellent agreement with 
experiment, especially for the lighter alkali chlorides (Li-K), there is a tendency for 
theory to underestimate D 0 and overestimate r e . Whereas the error for RbF is only 
0.017A , the error for RbCl is 0.043A at the CI(SD) level. For the heavier alkali 
halides where the bonding is predominantly ionic, it is tempting to use the nearly 

linear relationship between D e and r e to "correct” the theoretical D 0 values. Allow- 

0 SCF 

ing for a 0.03 eV change in D 0 per 0.01A change in r produces "corrected" D 0 

values of 4.37 eV for KC1 and RbCl and 4.54 eV for CsCl in excellent agreement 

with experiment. One produces similar results for KC1 and RbCl by correcting the 

D 0 results by 0.02 eV per 0.01A change in r, which is approximately the same 

correction that one would derive from the linear relationship among the alkali 

chlorides in Fig. 1. 

The "corrected" theoretical D 0 values provide strong support for the available 
experimental values. We prefer the thermochemical values of Brewer and Brackett 
[1] slightly over the flame photometry values of Bulewicz et al. [41], since our 
theoretical results suggest that the dissociation energy of RbCl is less than or 
equal to that for KC1. For CsCl we have not yet carried out the CI(SD) calcula- 
tion, but we estimate a value of 4.4 eV, significantly below the experimental values 
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[43—44 J of 4.58 eV. Whether this is entirely a result of basis set incompleteness or 
due partially to a tendency of the CI(SD) procedure to underestimate the pair-pair 
correlation, we hope to elucidate in future work [23]. 


D. Alkaline-earth chlorides 


The agreement between our CI(SD) D 0 values and experiment [54,60] is excel- 
lent if we take the lower mass spectrometric determination of 3.98±0.10 eV. [54] 
for BeCl. The very large positive correlation effect observed for BeCl arises from 
the large component of covalent character in the bonding. Note that BeCl has a 
r x /r e value of only 1.41 (see Table I). When there is a significant differential corre- 
lation effect, there is a tendency to underestimate it, so the mass spectrometric 
value of Hildenbrand [54] may be correct within their error bars. Thus, we recom- 
mend a value of 3.9±0.1 eV for BeCl. For MgCl, we also prefer the lower mass 
spectrometric determination of 3.26±0.13. The errors in the bond lengths for the 
heavier alkaline-earth chlorides are less than the corresponding alkali chlorides (see 
Table VII). Hence, the CI(SD) D 0 values for CaCl and SrCl are probably only 
slightly (0.05 eV) too low. Allowing for this the lower bound of 4.29 for SrCl deter- 
mined from chemiluminescent studies [60] would appear to be slightly high. 


E. Alkali oxides 


The alkali oxides are unique in that they undergo a change in ground state 
symmetry from 2 IT(LiO, NaO) to 2 E + (KO, RbO, CsO). This has been explained 
[15] in terms of the competing effects of quadrupole interactions (favoring 2 II) and 
Pauli repulsion (favoring 2 E + ). Both the theoretical [14-15] and experimental [12- 
13,64] evidence is that LiO and NaO have 2 II ground states, whereas RbO and 
CsO have 2 £ + groupd states. For KO there is conflicting experimental evidence; 
the ESR spectrum [13] suggests a 2 II ground state, whereas a magnetic analysis of 
the K + NO 2 system [12] suggests a 2 E + ground state. Previous theoretical results 



- 17 - 


for KO are conflicting as well, with the SCF calculations of So and Richards [14] 
giving a 2 E + ground state, and both the SCF and Cl calculations of Allison et. al. 
[15] giving a 2 II ground state. 

Our theoretical 2 II— 2 E + energy separations are compared with previous 
theoretical and experimental evidence in Table VIII. Throughout we have corre- 
lated 15 valence electrons (11 for LiO), which we believe is essential for a balanced 
treatment of the two states. For all of the alkali oxides we find that the 2 II— 2 E + 
energy separation is insensitive to electron correlation, because both states are 
equally well described by the HF configuration. This is also generally true of previ- 
ous theoretical results except for the KO results of Allison et al. [15] where they 
observe a 600 cm -1 increase in the separation. However, this increase is probably 
an artifact of correlating only the oxygen valence electrons without first localizing 
the orbitals. In KO there is a substantial mixing of the oxygen 2s and potassium 
3p orbitals, which is somewhat more pronounced in the Z E + state. The greater 
mixing in the z E' r state reduces its correlation energy with respect to the 2 II state 
producing the larger separation. When both the oxygen and metal electrons are 
correlated with the CI(SD) procedure, the computed energies are invariant to these 
arbitrary mixings. 

Apparently, basis set quality is a more important consideration than electron 
correlation for an accurate determination of the 2 II— 2 E + energy separation. This 
statement is substantiated by the results in Table DC, where we have studied basis 
set effects on the energy separation in KO at 4.2 Bohr. Starting with the large 
Slater basis set in Table IX, and selectively deleting functions, we find that delet- 
ing the diffuse functions has little effect on the energy, but deleting the 4f func- 
tions and particularly the 4f functions on oxygen, has a substantial (250 cm -1 ) 
effect on the separation. Apparently, f functions are more important for the 2 E + 
state. Similarly, deleting all d- and f-functions results in a very large (960 cm -1 ) 
effect on the 2 II— 2 E + separation, since the d functions are more important for the 
2 n state. Finally, to assess how close our Slater basis is to the HF limit, we have 
carried out numerical HF calculations on both the 2 IT and 2 E + states. The numeri- 
cal HF calculations (see Table DC) show that our Slater basis is within 0.018 and 
0.015 eV of the HF limit for the 2 £ + and 2 II states at 4.20 Bohr, respectively. We 
obtain a 2 II— 2 E + separation of -250 cm -1 in the HF limi t, which lies between our 
Slater result and that of So and Richards [14]. Hence, we agree with So and 
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Richards [14], and disagree with Allison et al. [15], who obtain a 2 IT ground state 
at the HF level, probably as a consequence of the relatively small basis set 
employed. Although our results are not definitive, they strongly support the 
assignment of a 2 E‘ r ground state for KO. However, the energy separation is prob- 
ably sufficiently small to allow the population of both states even at room tem- 
perature. 

The spectroscopic parameters (r e , w e ) for the 2 I1 and 2 £ + states of the alkali 
oxides are summarized in Table X. We observe almost identical bond contractions 
upon correlation in the 2 I1 states of the alkali oxides as for the 'E* states of the 
alkali fluorides. Hence, a comparison of the CI(SD) and experimental r e values for 
the alkali fluorides should provide a good measure of the remaining errors in the 
corresponding 2 I1 states of the alkali oxides. However, we observe a significantly 
larger CI(SD) bond contraction in the 2 E + states, probably because the hole in the 
2p^ orbital changes the character of the metal-oxygen repulsion. At the CI(SD) 
level we obtain nearly identical bond lengths for the 2 E + states of KO and RbO as 
for the 1 E + states of KF and RbF. This is significantly different from the assump- 

o 

tion of r e (MO) = r e (MF) +0.05A made by Herm and Herschbach [12] for purposes 
of estimating MO dissociation energies based on the corresponding alkali fluorides. 

Our theoretical dissociation energies for the alkali oxides (with respect to 
ground state atoms) are compared with experiment in Table XI. Our results gen- 
erally agree with the estimates of Herm and Herschbach [12], within their error 
bounds, illustrating again how well bond energies scale with bond length for these 
ionic systems. Our CI(SD) dissociation energies for LiO and NaO are significantly 
larger than the values derived from mass spectrometry [65-66]. In fact, they are 
larger by about the 2 II— 2 E + energy separation, which suggests the possibility that 
these molecules were prepared in the excited A 2 E + state. Our D 0 for LiO is also 
significantly larger than the theoretical value of Yoshimine [61], which was 
obtained by referencing directly to ground state atoms. This procedure underesti- 
mates the differential correlation energy thereby underestimating D 0 . We do, how- 
ever, agree with the theoretical estimate of O’Hare and Wahl [62] of 2.72±0.3 eV 
for NaO. Apparently, all of the dissociation energies reported by Allison et al. [15] 
are significantly too low as a result of basis set incompleteness accentuated, by 
their procedure of correlating only eight valence electrons. For KO we are in excel- 
lent agreement with the Ehlert [67] mass spectrometric determination, and agree 
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with the Gusarov and Gorokhov [68] value within their error bounds. 

For RbO( 2 £ + ) and CsO( 2 £ + ). there are to our knowledge no experimental 
values available for comparison, but we are in good agreement with the estimates 
of Herm and Herschbach [12]. Since our theoretical values tend to be slightly low 
for the heavier alkali atoms our recommended values are 2.9±0.1 eV for RbO and 
3.1±0.1 for CsO. We arrive at the value for CsO by adding to our SCF value a 
0.15 eV correction for electron correlation and a 0.05 eV correction for basis set 
incompleteness. 


F. Alkaline-earth oxides 


The 1 £ + ground states of the alkaline-earth oxides are not well represented by a 
single reference configuration reflecting the fact that the charge distribution is 
intermediate between M + 0* and . Therefore, we have applied our method 

to the a 3 IT states, (A 1 II for BeO) since they are well represented by the SCF refer- 
ence and the T g values are known experimentally [69]. For BeO we used the A 1 !! 
state since the T g for the a 3 II state is not accurately known. The spectroscopic 
parameters (r e , u; e ) are compared with the corresponding experimental values [69] 
in Table XII. The CI(SD) values are uniformly better than the SCF values, and 
show’ errors just slightly larger than for the fluorides. 

The CI(SD) D 0 values for the alkaline-earth oxides shown in Table XI should be 
the most consistent set of numbers available, and are thus capable of ruling out 
incorrect experimental values. Theoretical dissociation energies have been reported 
for both MgO [8] and CaO [ll] . The values reported herein are just slightly larger 
owing to the larger basis sets employed and the treatment of electron correlation. 

Our CI(SD) value of D 0 =4.69 eV for BeO is just slightly above the old mass 
spectrometric value of Chupka et al. [70] This value includes correlating the Be Is 
electrons, and produces a bond length for the A TI state that is 0.007A shorter 
than experiment. If the Be Is electrons are not correlated D o =4.60 eV and r e = 

c 

1.462 A in almost perfect agreement with experiment [69]. However, we believe 
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the higher value of D 0 =4.69 eV is preferred, and that the experimental values are 
just slightly too low. Note that we arrive at our values using a T e (A 1 II- 1 S + ) of 
9405.6 cm -1 , and correcting for the zero-point motion in the ground state using 
u,- e =1187.3 cm -1 . 

All of the available experimental values [2-7] for the 1 E + state of MgO are sub- 
stantially too large. The value recommended by Huber and Herzberg [69] is 3.5 eV, 
which is a reinterpretation of the Srivastava value [7] in Table XI, taking into 
account the presence of the low-lying a 3 II state. In contrast, the theoretical calcu- 
lations by Bauschlicher et al. obtain D o =2.65±0.16 eV. The CI(SD) value reported 
here of 2.75±0.1 eV is slightly larger as a result of the more extended basis set and 
Cl treatment. Note that this D 0 value is what would be expected based on the 
nearly linear relationship between bond distance and bond energy shown in Fig. 1. 
In determining D 0 calculations were performed for the a 3 II state, which was 
presumed to be 0.326 eV above the ground state [71]. Based on a comparison of 
our computed r e value for MgF^E*) with experiment, our calculations suggest 
that the somewhat uncertain r e value of 1.8 7 reported in Huber and Herzberg [69] 
is slightly long. 

For CaO, recent chemiluminescence studies have yielded different values for the 
dissociation energy. The study of Engelke, Sander and Zare (ESZ) [9] yields a 
lower bound of ^4.76±0.15 eV, while that of Irvin and Dagdigian (ID) [10] yields 
4.11 ±0.07 eV. The work of ID is in excellent agreement with the older mass spec- 
trometric values [3,7], but the ESZ result was recommended by Huber and 
Herzberg [69]. Recent SCF calculations by Bauschlicher and Partridge [11] gave 
strong support for the determination of ID. Our CI(SD) result of 4.14±0.1 is also 
in excellent agreement with the chemiluminescent value of ID. We thus recommend 
adopting this value for CaO, which again places the value on the linear plot of Fig. 
1. The higher value is now thought to be in error as a result of interference with 
CaCl. 

For the SrO molecule there is again an apparent discrepancy between the 
recent chemiluminescence studies of ESZ [9] that gave a lower bound of 
D o ^4.67±0.15 eV, and the older mass spectrometry and flame photometry values. 
Our value is in excellent agreement with the value recommended by Srivastava [7] 
in a review of the experimental literature prior to 1975. We expect our value is 
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slightly (~0.05 eV) too low judging from an overestimation of the bond length by 

c 

0.02A . Hence, we would recommend a value of 4.36±0.10 eV that rules out the 
determination by ESZ. The ESZ lower bound for SrO was based on the identifica- 
tion of a small spectral feature as the (18,0) band of the SrO system 

(see Fig 3 of [9]). However, owing to the strong interference from SrCl, all that 
can be safely determined is that the (11,0) band is populated. Since the (18,0) 
band is at 17,420 cm -1 and the (11,0) band is at 14,237 cm -1 , 3203 cm -1 should 
be removed from the lower bound, giving a revised estimate [72] of 
D o (SrO)^4.27±0.15 eV. 

For BaO we have not yet carried out the CI(SD) treatment, but have deter- 
mined an SCF value of 5.61 eV. Since we anticipate that D 0 will increase by at 
least 0.1 eV with the CI(SD) treatment as well as with further basis set saturation, 
we recommend a value of 5.75±0.10 eV. Although the error bars given for BaO 
may be somewhat optimistic, our recommended value is in good agreement with 
the lower bound of ^5.79±0.15 eV given by ESZ [9] determined from the chemi- 
luminescence spectrum of Ba + C10 2 . Note also that their determination for BaO 
seems much more definitive since there is no interference from BaCl (see Fig. 2 of 

l®])- 


G. Alkali sulfides 


We are unaware of any experimental work on the alkali sulfides (MS) which 
are difficult to observe as a result of their strong tendency to form M 2 S. In Table 
XIII we have summarized our theoretical MS(X 2 II) D 0 values and X 2 IT-A 2 E + 
energy separations for the alkali sulfides LiS-RbS. The ground states for all of 
these systems are definitively determined to be X 2 II. Although calculations have 
not yet been done for CsS, we predict a 2 II ground state as well. The 2 II states 
would appear to be more stable with respect to the 2 E + states in the sulfides than 
the oxides, because the considerably longer bond lengths (r e (MS)«1.22r e (MO), 
M=Na, K, Rb) reduce the importance of the Pauli repulsion terms. 

The spectroscopic parameters (r e , u/J for the X 2 II and A 2 E + states of the alkali 
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sulfides (Li-Rb) are given in Table XIV. For the X 2 II state we observe a nearly 
identical bond contraction upon correlation as observed for the 1 E + states of the 
corresponding alkali chlorides. For the A 2 E^ states we observe a larger Cl bond 
contraction than in the 2 IT states, in analogy with the alkali oxides. 

Since the X 2 II states of the alkali sulfides have considerable analogy with the 
X J E + states of the corresponding alkali chlorides, we expect the same deficiencies 
in the CI(SD) D 0 values and bond lengths of the heavier alkali sulfides (KS, RbS) 
as for the chlorides. Thus, in determining our recommended D 0 values in Table 
XIX, we have added 0.06 and 0.08 eV, respectively, based on the known deficien- 
cies in our D 0 values of the corresponding alkali chlorides. Note that the alkali sul- 
fides obey the nearly linear relationship between dissociation energy and bond 
length illustrated in Fig 1, even though by the r x /r e criterion they are far less ionic 
than the corresponding chlorides. 


H. Alkaline-earth sulfides 


For the alkaline-earth sulfides (BeS-CaS) we have encountered the same prob- 
lem as for the oxides, namely that the X 1 E + ground states are not well described 
by a single reference configuration, and therefore not amenable to our model. 
Therefore, we have considered instead the A 1 II and a 3 II states, which are both 
equally well described by a single reference configuration. Our theoretical spectros- 
copic parameters (r e , u> e , D e ) for the A 1 !! states of BeS, MgS and CaS as well as 
the a 3 II-A 1 n energy separations are summarized in Table XV. The corresponding 
spectroscopic parameters (r e , u/J for the a 3 II state (not reported) are essentially 
identical to those for the A 1 !! state. The theoretical a 3 II-A 1 n energy separations 
for the sulfides decrease down the column as the K(<t,tt) exchange integral 
decreases with increasing internuclear distance since the a open-shell orbital is 
localized on the metal, and the ir electron is primarily a 3p7r sulfur orbital. The 
singlet-triplet separations are comparable or slightly smaller for BeS and MgS than 
for the corresponding oxides [73-76]. 

To convert our dissociation energies for the A 1 IT states in Table XV to ground 
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state D 0 . we must add T e (A 1 II) and subtract the X J E + zero-point energy. Unfor- 
tunately, T (.A^II) is known only for BeS, and even here the vibrational numbering 
in the A 1 II state is uncertain by one unit [77]. Using T e (A 1 II) = 7842.9 cm -1 and 
u.’ e (X 1 E‘ r )=997.94. our D 0 (X 1 H + ) for BeS is 3.20 eV. The theoretical value is con- 
sistent with the rather uncertain value of 3.47±0.65 eV reported in the JANAF 
tables [78]. A linear Birge-Sponer extrapolation of the ground state with correc- 
tions for ionic character [79] and excited state products yields a value of 3.30 eV 
[78). 

Our theoretical dissociation energy (D e ) for the A 1 II state of MgS is 1.70 eV. If 
we assume in analogy with MgO that the A 1 II state is low-lying (assuming 
T e =3000 cm -1 following JANAF, we obtain D o (X 1 E + )=2.04 eV, which is consider- 
ably lower than the rather uncertain JANAF value [78] of 2.86±0.69 eV. Hence, it 
is likely that like MgO, the JANAF value for MgS is significantly too large. 

Our CI(SD) dissociation energy for the A 1 !! state of CaS is 2.41 eV. Adopting 
the rather uncertain JANAF estimate of 7200 cm -1 for T e (A*n), produces a 
D 0 (X 1 £' 1 ') of 3.27 eV in excellent agreement with the mass spectrometric studies of 
Colin et al. (D o =3.20±0.20 eV) [80], but less than that of Marquart and Berkowitz 
(3.45) [81]. This suggests that the a 3 II and A 1 !! excitation energies for the 
alkaline-earth sulfides are reasonably similar to the corresponding oxides. Theoreti- 
cal work is in progress to accurately determine these excitation energies [26], 


I. Alkali hydroxides 


Our theoretical dissociation energies for the alkali and alkaline-earth hydrox- 
ides are compared to experiment [82-97] in Table XVI. The most striking feature 
of the alkali hydroxides (especially the heavier ones) is their similarity to the alkali 
fluorides. The dissociation energies of the hydroxides are less than the fluorides by 
slightly less than the electron affinity difference (EA(F)-EA(OH)=1.57 eV) [21,98]. 

Comparison of the alkali hydroxides with the alkali oxides is also enlightening. 
Here the electron affinity difference (EA(OH)-EA(O)~0.37) is only about half of 
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the difference in dissociation energies (D 0 (MOH)-D 0 (MO)). The increased stability 
of the hydroxides relative to the oxides is a result of the ability of the hydrogen 
atom to pull charge density out of the metal oxygen bond, thereby reducing the 
bond length and giving rise to greater electrostatic interaction [99]. The increased 
stability of the alkali hydroxides with respect to the oxides is also evident from the 
plot of D e versus r e in Fig. 1. 

The theoretical CI(SD) D 0 values for the alkali hydroxides are in reasonably 
good agreement with the two sets of experimental values derived from atomic 
absorption spectroscopy in flames. The theoretical values are probably more accu- 
rate than the flame values, in part because the experimental values are based on 
equilibrium constants, which in turn require a knowledge of the rotational and 
vibrational partition functions of MOH. Our theoretical values for the hydroxides 
should be almost as accurate as for the fluorides, and show an amazing parallel. 
There is a tendency for theory to overestimate the bond length slightly for the 
heavier alkali hydroxides. Hence, based on the nearly linear relationship between 
bond length and D e , our recommended values in Table XIX add 0.02, 0.04 and 
0.06 to our CI(SD) D„ for KOH, RbOH and CsOH, respectively. These changes 
improve agreement with experiment, but our recommended value for CsOH of 
3.77±0.10 eV is still considerably less than the experimental values [82,84]. Note 
that the difference in dissociation energy of CsOH and RbOH obtained from the 
two flame experiments was quite different (0.13 eV versus 0.34 eV). We obtain 
about 0.19 eV for this difference, which is essentially the same as the difference in 
dissociation energy of CsF and RbF. 

Our theoretical r e and u; e values for the alkali hydroxides are compared with 
experiment [100-109] in Table XVII. These results were obtained assuming a fixed 

o 

OH bond length of 0.947A . If we correct our CI(SD) r e values for the hydroxides 
by the known deficiencies in the CI(SD) values of the corresponding fluorides, then 
we are in satisfactory agreement with the rather uncertain experimental values, 

o 

except for NaOH. For NaOH we estimate a Na-0 bond length of 1.936A , which is 

o 

significantly shorter than the experimental estimate [101] of 1.95A . Our theoreti- 
cal ui e values for the alkali hydroxides tend to be slightly larger than the 
corresponding fluorides, arising primarily from the difference in mass, i.e. (Mq + 
Mjj) < Mp. The theoretical uj e values are probably the most accurate available. 
For CsOH we agree with the rather uncertain gas phase value [106]. However, our 
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results imply that all of the u e values derived from matrix isolation studies 
[102,104,107] are systematically low. 


J. Alkaline-earth hydroxides 


The similarity of the alkaline-earth hydroxides and fluorides is demonstrated 
by the recent high-resolution laser excitation spectra of CaOH [108] and SrOH 
[109]. The lowest electronic transitions for both the alkali fluorides and hydroxides 
involves the promotion of a nonbonding electron localized primarily in the metal 
ns<r orbital to a metal np?r or np<7 orbital. Since these transitions are localized pri- 
marily on the metal, the potential curves for each state are similar, giving a spec- 
trum consisting of very strong, badly overlapped Av=0 sequences. 

The theoretical dissociation energies of the alkaline-earth hydroxides are smaller 
than the corresponding fluorides by less than the difference of 1.57 in the electron 
affinities of F and OH. Again this demonstrates the increased stability of the 
hydroxides that arises from the ability of hydrogen to pull charge, out of the 
metal-oxygen bond. Experimental studies of the alkaline-earth hydroxides in flames 
are complicated by the presence of the dihydroxide. In general, the flame-spectral 
data appear to give dissociation energies that are too high. The experimental 
values for BeOH are particularly uncertain. Our D 0 value of 4.70 eV for BeOH is 
within the error bounds of the value of 4.94±0.43 reported in the JANAF Tables 
[87]. which is the average of two values ascribed to Inami and Ju [86]. The consid- 
erably larger D 0 value of Ko et al. [85] for BeOH was discarded based on a com- 
parison of the trends in the alkaline-earth hydroxides and halides. Note that our 
value of D o =4.70 eV was determined for a linear geometry and includes correlating 
the Be Is electrons. However, calculations with a smaller gaussian basis set lead to 
a bent (0=33°) equilibrium structure about 0.01 eV more stable than the linear 
one. 

For MgOH, the two older experimental values of D o =2.43±0.22 eV derived by 
Bulewicz and Sugden [89] from flame spectra studies, and D o =4.16±0.13 eV 
derived by Cotton and Jenkins [88] using atomic absorption spectroscopy appear 
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to be in error. The JANAF Tables [87] have adopted an intermediate value of 
3.59±0.22 eV, which is based primarily on the trends between the hydroxides and 
halides. There is also a more recent mass spectrometric determination [90] of 
3.19±0.22 eV that is more reliable. Our theoretical value of D 0 =3.3l eV is some- 
what larger, but agrees with the mass spectrometric determination within its error 
bars. 

For the heavier alkaline-earth hydroxides, there is a larger body of reliable 
experimental D 0 values. Our theoretical D 0 for BaOH appear to be systematically 
low by ~0.2 eV as a result of using a relativistic effective core potential. Hence for 
BaOH we are only able to recommend a value of D o =4.6±0.2 eV, which includes 
all of the experimental values [90,93,97] except that of Cotton and Jenkins [92], 
which appears to be too high as are their values for CaOH and SrOH. 

For CaOH and SrOH our CI(SD) D 0 values agree with the mass spectrometric 
determinations [90] within their error bounds, but strongly suggest that the values 
are systematically 0.10-0.15 eV low for the hydroxides of Mg, Ca and Sr. The 
lower flame values of Gurvich et al. [96] for SrOH and Ryabova et al. [93] for 
CaOH agree quite well with our theoretical values. Overall, our theoretical values 
should be the most reliable, and may be capable of shedding light on the relative 
equilibrium constants for the formation of the mono- and di-hydroxides. 

The theoretical spectroscopic constants r e and w e for the alkaline-earth hydrox- 
ides are summarized in Table XVII. These were determined using a Fixed OH bond 

o 

length of 0.947A . For BeOH we obtain a non-linear equilibrium structure with an 
equilibrium bond angle of about 33°. The heavier alkaline-earths all favor linear 
structures, but are exceptionally flat in the bending potential. The u; e values for 
the alkaline-earth hydroxies are about 10% greater than the corresponding 
fluorides. Our w e values for CaOH and SrOH are in very good agreement with the 
gas phase values [108,109]. Little is known about the equilibrium bond lengths for 
the alkaline-earth hydroxides. Our CI(SD) values for CaOH and SrOH axe slightly 
larger than the gas phase values [108,109] as expected. However, our r e values for 
BeOH and MgOH should be quite accurate (see corresponding fluorides). 
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K. Alkaline-earth isocyanides 


We summarize here the theoretical calculations [l 10] that were undertaken to 
determine the lowest energy structures and dissociation energies for the 2 A' 
ground-state surfaces of MCN (M=Be, Mg, Ca and Ba). There has been consider- 
able interest in the corresponding alkali cyanide molecules recently with the 
discovery that both NaCN and KCN have non-linear T-shape equilibrium struc- 
tures [111-115]. The dominant features in the spectra of the alkaline-earth 
(iso)cyanides are broad "quasicontinua", that occur at nearly the same 
wavelengths as the spectra of the homologous alkaline-earth monohalides [116]. 
These molecules are also expected to be ionic based on the successive change in 
electronegativities C1<CN<F. Note, however, that the electron affinity [117] of 
CN, 3.82±0.02, is considerably greater than even that of fluorine (3.399 eV). 

For all of the alkaline-earths studied, the linear isocyanide structure was found 
to be the most stable. At the CI(SD) level, the cyanide structure was found to lie 
above the isocyanide structure by 0.26 eV for Be, 0.13 eV for Mg and 0.20 eV for 
Ca. One interesting change in going from BeCN to BaCN is the loss of the inter- 
conversion barrier between the isocyanide and cyanide structures. On the basis of 
several cuts through the potential surface, we found an interconversion barrier for 
BeCN and MgCN, but none for the heavier alkaline-earths, where the bonding 
becomes increasingly ionic. 

Our theoretical dissociation energies for the cyanide, isocyanide and bond-mid- 
point (metal atom approaching CN bond midpoint) structures are compared with 
experiment in Table XVIII. The experimental dissociation energies [118] for the 
monocyanides have been determined using electrothermal atomic absorption spec- 
trometry. These values contain considerable uncertainty, since little was known 
about the number of low-lying electronic states, geometry or force constants of 
these molecules. The experimental values lie consistently above our values, particu- 
larly for MgCN w’hich we feel must contain some systematic error. On the other 
hand, our value for the BaCN structures are probably 0.2 eV too low as a result of 
using a relativistic effective core potential treatment of the core. 


Our dissociation energies for the isocyanides are very similar to the 
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corresponding chlorides. This correlation provides additional support for our 
theoretical value for MgCN. Apparently, the delocalization of charge in CN makes 
the D e of the isocyanides more comparable to the chlorides, even though the M-NC 
bond lengths are intermediate between the corresponding alkaline-earth fluorides 
and chlorides. 
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I. Conclusions 


The theoretical dissociation energies presented here for most of the alkali and 
alkaline-earth fluorides, chlorides, oxides, sulfides, hydroxides and isocyanides are 
sufficiently accurate to rule out disparate experimental values. Overall, we find 
that the thermochemical and mass spectrometric determinations of D 0 are quite 
accurate. The D 0 values determined by flame photometry are often accurate, but 
tend to be systematically high (especially for the fluorides). Finally, we find that 
D 0 values determined from chemiluminescent studies, although potentially very 
accurate, are often incorrect. 

An advantage of treating several classes of systems in a systematic way, is that 
trends in the dissociation energies emerge more clearly. For example, the MF and 
MOH bond lengths are comparable, and the difference in dissociation energies is 
slightly less than the electron affinity difference of 1.57 eV. However, the iso- 
cyanides have distinctly smaller dissociation energies than the fluorides even 
though the electron affinity of CN is substantially greater than that of F. The 
difference arises because the charge is spread out in CN, but localized on oxygen in 
OH. 

The theoretical D 0 values in conjunction with the available experimental values 
allow us to recommend reliable D 0 values for all systems considered. These recom- 
mended values in Table XIX are thought to be accurate to 0.1 eV. The model we 
have applied seems to be satisfactory even for systems (e.g. BeCl) that are not 
fully ionic as along as we account for differential correlation effects with the 
CI(SD) procedure. We are willing to calculate the dissociation energy of any other 
diatomic or triatomic system that fits our model, if it is of sufficient scientific 
interest. 
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Table I. Ratio r x /r e for selected alkali and alkaline earths 


Molecule 

r e (A ) a 

r,(A ) b 

Ratio 

LiF( ! E + ) 

1.564 

7.232 

4.62 

LiCl( 1 E + ) 

2.021 

8.113 

4.01 

LiO( 2 n) 

1.695 

3.666 

2.16 

LiS( 2 n) 

(2.163) 

4.347 

2.01 

NaF( 1 E + ) 

1.926 

8.280 

4.30 

NaCl( 1 E + ) 

2.361 

9.455 

4.01 

NaO( 2 II) 

(2.040) 

3.917 

1.92 

NaS( 2 n) 

(2.487) 

4.705 

1.89 

kf(*e + ) 

2.171 

15.319 

7.06 

KC1( 1 E + ) 

2.667 

19.889 

7.46 

KO( 2 E + ) 

(2.187) 

5.005 

2.29 

KS( 2 n) 

(2.843) 

6.367 

2.24 

RbF( 1 E + ) 

2.270 

18.533 

8.16 

RbCl^E 4 ) 

2.787 

25.668 

9.21 

RbO( 2 E + ) 

(2.288) 

5.306 

2.32 

RbS( 2 n) 

(2.983) 

6.861 

2.30 

BeF( z E + ) 

1.361 

2.432 

1.79 

BeCl( 2 E + ) 

1.797 

2.524 

1.41 

BeO( 1 E + ) 

1.331 

1.832 

1.38 

BeS( l E + ) 

1.742 

1.988 

1.14 

MgF( 2 E + ) 

1.750 

3.392 

1.94 

MgCl( 2 E + ) 

2.199 

3.574 

1.63 

MgO( J E + ) 

1.749 

2.329 

1.33 

MgS( 1 S + ) 

2.143 

2.587 

1.21 

CaF( 2 E + ) 

1.967 

5.310 

2.70 

CaCl( 2 E + ) 

2.439 

5.769 

2.37 

CaO( 1 E + ) 

1.822 

3.097 

1.70 

CaS( 1 E + ) 

2.318 

3.570 

1.54 

SrF( 2 E + ) 

2.075 

6.280 

3.03 

SrCl( 2 E + ) 

2.576 

6.933 

2.69 

SrO^E*) 

1.920 

3.404 

1.77 

SrS( 1 E + ) 

2.441 

3.984 

1.63 
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a Bond distances in parentheses are theoretical values for cases where experimen- 
tal values are unavailable. 

^Hypothetical crossing point where a purely ionic potential curve (1/R) crosses 
the ground state asymptote. 
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Table II. Summary of Slater and gaussian basis sets 


A. Slater Basis sets 


Atom 

Basis size 

Basis set description a 

O 

7s5p3dlf 

dementi^ (6s4p) O and 0“ basis. Polariza- 
tion functions optimized for O and 0" 

F 

6s5p4d2f 

dementi (6s4p) F and F~ basis. Polariza- 

d 

tion functions optimized for F and F~ 

S 

9s7p3d2f 

dementi (7s6p) S“ basis. Even tempered d 
and f functions optimized at the CI(SD) 
level. 

Cl 

9s7p3d2f 

dementi (7s6p) Cl“ basis. Even tempered 
d and f functions optimized at the CI(SD) 
level. 

Li 

6s6p4d2f 

Konowalow e (5s5p3d) basis 

Li(ls) 

6s8p6d2f 

Adds compact 2p (4.0, 8.0) and 3d (5.0, 
8.0) 

Na 

8s6p4d2f 

dementi (7s, 3p) basis 

K 

10s8p4d2f 

dementi (lls6p) basis with two most dif- 
fuse s functions deleted 

Rb 

Ils9p6d2f 

dementi (Ils7p3d) Rb basis minus most 
diffuse s 

Cs 

14sllp9d2f 

McLean^ (I0s8p4d) Cs~*~ basis reoptimized 

Be 

6s4p3dlf 

Liu^ Be 2 basis 

Be(ls) 

6s6p5d2f 

Adds 2p(6.0, 3.5), 3d(8.0, 4.0) and 4f(2.8) 

Mg 

8s6p4d2f 

8s6p4dlf basis used for MgC> 

Ca(3d) 

I0s8p3d2f 

Ca basis used for CaO* 3d(0.8007, 2.122, 
4.77) 

Ca(4d) 

10s8p4d2f 

Even tempered d basis 3d(0. 75, 1.5, 3.0, 
6.0) 
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Sr 13sl0p7d3f dementi (Ils7p3d) Sr + basis reoptimized. 

Valence exponents optimized for the Sr P 
and states. 

Ba 15sl3p9d3f McLean (12s8p4d) Ba + basis reoptimized. 

3 

Valence exponents optimized for the Ba P 
and states. 

B. Gaussian basis sets 

O (Ils7p3dlf/6s4p3dlf) van Duijneveldt^ lls6p. 

H (8s4p/5s3p) van Duijneveldt 6s. 

C (Ils7p2d/6s4p2d) van Duijneveldt lls6p. 

(Ils7p3dlf/6s4p3dlf) Expanded d and f basis. 

N (Ils7p2d/6s4p2d) Vein Duijneveldt lls6p. 

(Ils7p3dlf/6s4p3dlf) Expanded d and f basis. 

Li (Ils7p7d4f/7s7p7d3f) van Duijneveldt 11s. 

Be (Ils5p2d/7s3p2d) van Duijneveldt 11s. 

(Ils6p7d4f/7s6p7d3f) Tight p and expanded d and f basis. 

Na (12s9p6d3f/7s6p5d2f) McLean and Chandler^ 12s9p. 

Mg (12s9p5d/6s6p4d) McLean and Chandler 12s9p. 

(12s9p6d4f/6s6p5d3f) Expanded d and f basis. 

K (15sl3p6d4f/9s9p5d3f) Wachters* 14s9p. 

Ca (15sl3p5d/9s9p3d) Wachters 14s9p. 

(15sl3p7d2f/9s9p6d2f) Expanded d and f basis. 

Rb (18sl4plld4f/12sl0p5d3f)Huzinaga m 17sllp6d. 

Sr (17sl4plld4f/lls9p5d2f) Huzinaga 17sllp6d. 

Cs (14s9p6d4f/7s7p5d3f) RECP of Christiansen and Laskowski 11 , 

modified valence basis. 

Ba (8s9p6d/6s7p4d) RECP of Mascarello and Jaffe°. 

(8s9p6d4f/7s7p4d3f) Recontracted s basis and added f functions. 

a A full description of the basis sets is given the the complete works, [23-28]. 
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^E. Clementi and C. Roeti, At. Data Nucl. Data Table 14, 177 (1974). 

c C. W. Bauschlicher, unpublished. 

^B. Liu, private communication. 

e D. D. Konowalow, and M. L. Olson, J. Chem. Phys., 71, 450 (1980). 

*A. D. McLean and R. S. McLean, IBM Research Report RJ 3187 (1981). 

g B. H. Lengsfield, A. D. McLean, M. Yoshimime and B. Liu, J. Chem. Phys., 79, 
1891 (1983). 

h Ref. [8]. 

'Ref. [11]. 

■^F. B. van Duijneveldt, IBM Research Report RJ 945 (1971). 

^A. D. McLean and G. S. Chandler, J. Chem. Phys., 72, 5639 (1980). 

! A. j. H. Wachters, J. Chem. Phys., 52, 1033 (1970). 

m S. Huzinaga, J. Chem. Phys., 66, 4245 (1977). 

n B. C. Laskowski, S. P. Walch and P. A. Christiansen, J. Chem. Phys., 78 
(1983) 6824, modified valence basis tabulated in, S. R. Langhoff, C. W. 
Bauschlicher and H. Partridge, J. Phys. B, J. Phys. B: At. Mol. Phys. 18, 13 
(1985). 

°C. W, Bauschlicher, R. L. Jaffe, S. R. Langhoff, F. G. Mascarello, and H. Par- 
tridge, J. Phys. B, in press; “Oscillator strengths of some Ba lines; a treatment 
including core-valence correlation and relativistic effects”. 



- 2 - 


Table III. Correlation effects on D e and r in the X 2 S + 

e e 

ground state of NaF 


Description 

D e (eV) 

r e (A°) 

HF 

4.957 

1.924 

CI(SD) 8 electrons on F a 

4.876 b 

1.944 

CI(SD) 8 electrons on Na a 

4.956 

1.934 

CI(SD) 8 electrons on F 
+pair-pair terms 3 

5.034 

1.912 

CI(SD) 16 electrons 

5.023 

1.922 


a Orbitals were first localized using corresponding orbitals [30] by maximizing 
the overlap with the Na+ orbitals. 

k D e is 4.83 eV if the orbitals are not localized. 
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Table IV. Dissociation energies (D 0 ) of the alkali 
and alkaline-earth fluorides 


Molecule 

Theoretical D 0 (eV) a 

Experimental D 0 (eV) 

LiF 

6.05(6.06) 

5.90±0.35 b , 5.96 c 

NaF 

4.92(4.99) 

5.25±0.30 b , 4.94 c , ^5.33±0.13 d 

KF 

4.96(5.07) 

5.07±0.35 b , 5.10 c 

RbF 

4.90(5.04) 

5.20±0.35 b , 5.04 c 

CsF 

5.06( ) 

5.33±0.35 b , 5.19 c , 5.27±0.06 e , 
5.32±0.08 f 

BeF 

5.88(5.94) 

5.85±0.10 g , 6.25 h 

MgF 

4.64(4.66) 

4.79±0.05 1 , 4.62±0.10j 
5.4±0.2 k , S.SSiO.SO 1 

CaF 

5.43(5.51) 

SrF 

5.47(5.58) 

5.62±0.07\ 5.43±0.10^, 5.72±0.30 


a SCF results given first with CI(SD) results in parentheses. Zero-point correc- 
tion is one-half corresponding u/ e . 

^Bulewicz, Phillips and Sugden(196l) [41], flame photometry. 

c Brewer and Brackett (1961) [l], thermochemical(298K) . 

^Ham(1974) [42], chemiluminescent studies. 

e Berkowitz(197l) [43], photoionization. 

f 

Parks and Wexler(l984) [44], collision induced dissociation. 
s Hildenbrand and Murad(1966) [46], mass spectrometry. 

^Farber and Srivastava(1974) [47], mass spectrometry. 

1 Hildenbrand(1968) [49], mass spectrometry at 298K. 

J Ehlert, Blue, Green and Margrave(1964) [48], mass spectrometry at 298K. 

i. 

Blue, Green, Bautista and Margrave(1963) [50], thermochemical at 298K. 
^Ryabova and Gurvich(1964) [52], equilibria in flames. 
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Table V. Theoretical spectroscopic parameters (r e , u> e ) for the 
ground states of the alkali and alkaline-earth fluorides 


Molecule 


r e (A ) b 



u/ e (cm 

,) b 


SCF 

CI(SD) 

Expt c 

SCF 

CI(SD) 

Expt c 

LiF 

1.554 

1.557 

1.564 

926 

937 

910 

NaF 

1.924 

1.922 

1.926 

534 

531 

536 

KF 

2.205 

2.185 

2.171 

402 

437 

428 

RbF 

2.311 

2.287 

2.270 

367 

377 

376 

CsF 

2.400 


2.345 

327 

352 


BeF 

1.353 

1.358 

1.361 

1220 

1224 

1247.36 

MgF 

1.740 

1.746 

1.750 

731 

729 

711.69 

CaF 

1.980 

1.966 

1.967 

611 

595 

581.1 

SrF 

2.104 

2.086 

2.075 

478 

524 

502.4 


a The ground state is 1 £ + for the alkali fluorides and 2 E + for the alkaline-earth 
fluorides. 

^Spectroscopic parameters are based on a parabolic fit to the lowest energy 
points on a 0.1 Bohr grid. 

c Values taken from Huber and Herzberg [69j. 
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Table VI. Dissociation energies (D 0 ) of the alkali and 
alkaline- earth chlorides 


Molecule 

Theoretical D 0 (eV) a 

Experimental D 0 (eV) 

LiCl 

4.79(4.89) 

4.79±0.13 b , 4.85 c 

NaCl 

4.08(4.21) 

4.23±0.09 b , 4.23 c 

KC1 

4.15(4.31) 

4.32±0.09 b , 4.39 c 

RbCl 

4.10(4.28) 

4.40±0.09 b , 4.37 c 

CsCl 

4.19( ) 

4.65±0.13 b , 4.61 c , 4.58±0.07 d 
4.58±0.08 e 

Bed 

3.57(3.87) 

3.98±0.10 f , 4.51 g 

MgCl 

3.08(3.24) 

3.26±0.13 h , 3.36±0.02 1 

CaCl 

3.92(4.10) 

4.09±0.13 h , 4.15±0.13^, 4.08±0.07 
4.16±0.13 h , 4.15±0.08 k , ^4.29* 

SrCl 

4.04(4.17) 


SCF results given first with CI(SD) results in parentheses. 
^Bulewicz, Phillips and Sugden(1961) [41], flame photometry. 
c Brewer and Brackett(1961) [l], thermochemical(298K). 
^Berkowitz(1969) [43], photoionization mass spectrometry. 
Parks and Wexler(1984) [44], collision-induced dissociation. 

f 

Hildenbrand and Theard(1969) [54], mass spectrometry. 
s Farber and Srivastava(1973) [55], mass spectrometry. 
^Hildenbrand(l970) [56], mass spectrometry, 
barber and Srivastava(1976) [57], mass spectrometry. 
J Zmbov(l969) [58], mass spectrometry. 

jj 

Gurvich, Ryabova and Khitrov(1974) [59], flame photometry. 
^ Jonah and Zare(1971) [60], chemiluminescent. 


* 
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Table VII. Theoretical spectroscopic parameters (r e , u> e ) for the 
ground states of the alkali and alkaline- earth chlorides 


Molecule 3 


r e (A ) b 



w e (cm 1 ) b 



SCF 

CI(SD) 

Expt 

SCF 

CI(SD) 

Expt c 

LiCl 

2.036 

2.018 

2.021 

675 

668 

643 

NaCl 

2.390 

2.366 

2.361 

362 

354 

366 

KC1 

2.740 

2.698 

2.667 

258 

273 

281 

RbCl 

2.877 

2.830 

2.787 

223 

235 

228 

CsCl 

3.024 


2.906 

196 


214.17 

BeCl 

1.812 

1.797 

1.797 

865 

853 

846.7 

MgCl 

2.212 

2.199 

2.199 

478 

469 

462.12 

CaCl 

2.502 

2.458 

2.439 

355 

322 

367.53 

SrCl 

2.645 

2.607 

2.576 

287 

304 

302 


a The ground state is for the alkali chlorides and 2 E + for the alkaline-earth 
chlorides. 

^Spectroscopic parameters are based on a parabolic fit to the lowest energy 
points on a 0.1 Bohr grid. 

c Values taken from Huber and Herzberg [69]. For SrCl we used the value 
ascribed to T. Torring in Table A7.2 of the Ph. D thesis of P. Bernath(MIT). 
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Table VIII. 2 II— 2 E + excitation energies of the alkali oxides 


Molecule 

2 n 


This work a 

LiO 

2359(2391) 

NaO 

1429(1701) 

KO 

-205(-240) 

RbO 

-516(-650) 

CsO 

-798 ( ) 


2 n- 2 S + excitation energies (cm ] ) 


Other Theory 

2894(2634) b , 2342(2330) c 
2088(2177) b , 1236°, 1785 e 
233(831) b -347 d 
-138(-114) , -606 d 
-735(-846) b , -497(-726) f 


Expt 

>0 g 

>0 h 

<(P 

<o h -j 


HF results given first with the Cl results in parentheses, 
k (SCF /8-electron Cl) results of Allison et al. [15]. 
c (SCF/CI) results of Yoshimine [61]. 

^ SCF results of So and Richards [14]. 
e SCF results of O’Hare and Wahl [62]. 

^ (SCF/CI) results of Laskowski et al. [63]. 
g Freund et al. [64]. 
k Herm and Herschbach [12]. 

1 Evidence has been presented separately for both a 2 E + [12], and a 2 II [13] 
ground state. 

J Lindsay et al. [13]. 



Table IX. Basis set study of KO at the HF level a 


Basis set 

Energy (eV)^ 

Description 

2 n 

2 U + 

NHF 

-.015 

-.018 

[I8sl4p8d4f] b 

0.0 

0.0 

-diffuse functions(df) 

0.0035 

0.0030 

-df and 0(4f) 

0.0042 

0.029 

-all f-functions 

0.020 

0.052 

-all f-functions 

0.359 

0.320 

and all d-functions on K 
-all d- and f-functions 

0.608 

0.489 


a All calculations done at 4.20 Bohr. 

k All energies are given relative to the [I8sl4p8d4f] Slater calculation. The basis 
set consists of the potassium and oxygen basis sets in Table II plus a set of dif- 
fuse s, p, d and f functions on oxygen. The numerical Hartree Fock (NHF) 
energies at 4.2 Bohr are -674.01582 a.u. for the 2 E + state and -674.01249 a.u. 
for the II state. 
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Table X. Theoretical spectroscopic parameters (r e , w e ) for the 
2 n and 2 E + states of the alkali oxides® 

Molecule r e (A ^(cm -1 )^ 



2 n 

2 S + 

2 n 

2 S + 

LiO 

1.674(1.677) 

1.579(1.574) 

816(805) 

767(757) 

NaO 

2.046(2.040) 

1.953(1.943) 

519(511) 

553(543) 

KO 

2.351(2.329) 

2.225(2.187) 

409(396) 

417(439) 

RbO 

2.466(2.443) 

2.334(2.288) 

328(349) 

366(380) 

CsO 

2.562( ) 

2.426( ) 

327( ) 

331 ( ) 


Spectroscopic parameters are based on a parabolic fit to the lowest energy 
points on a 0.1 Bohr grid. 

^The SCF values are given first with the CI(SD) values given in parentheses. 
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Table XI. Dissociation energies (D 0 ) of the alkali and 

alkaline- earth oxides 


Molecule Theoretical D 0 (eV) a ' c Experimental D 0 (eV) 


LiO( 2 n) 

3.79(3.78) 

NaO( 2 II) 

2.75(2.83) 

KO( 2 E + ) 

2.76(2.86) 

RbO( 2 E + ) 

2.74(2.88) 

CsO( 2 E + ) 

2.90(- — ) 

BeO 

4.61(4.69) 

MgO 

2.71(2.75) 

CaO 

4.07(4.14) 

SrO 

4.23(4.32) 

BaO 

5.61(- — ) 


3.56±0.17 d , 3.49±0.06 e 

2.90±0.13 d , 2.61±0.20 f 

3.08±0.13 d , 2.86±0.13 s , 3.08±0.26 h 

2.95±0.13 d 

3.03±0.13 d 

4.6±0.1 1 , 4.53±0.13^ 

3.71±0.13^, 3.73±0.22 k 

3.95±0.09^, ^4.76±0.15 1 , 4.11±0.07 m 

4.03±0.22 k 

4.27±0.09^, ^4.27±0.15 1,n 
5.62±0.04^, £5.79±0.15 J 


2L 

SCF values are given first with CI(SD) values given in parentheses. 

^For the alkali oxides we report D c values for the designated ground state with 
respect to ground state atoms even though the 2 E + states dissociate to M( 2 P) 
+ 0( 3 P). Zero-point corrections are made using the calculated w e values. 

For the alkaline-earth oxides, D 0 values are reported for the X E state with 
respect to ground state atoms. Zero-point corrections are made using the exper- 
imental u} e values. 

^Estimate of Herm and Herschbach, [12] based on the alkali halide molecules. 
e Hildenbrand(1972) [65], electron-impact mass spectrometry. 

f 

Hildenbrand and Murad(1970) [66], mass spectrometry. 

®Ehlert(1977) [67], mass spectrometry at 298K. 

n Gusarov and Gorokhov(1971) [68], third-law treatment of mass-spectrometric 
effusion data. 

'Chupka, Berkowitz and Giese(1959) [70], mass spectrometry. 

J Srivastava(1976) [7], recommended values- review of experimental literature 
through 1975. See also papers by Farber and Srivastava [6,47,55,57]. 

k 

Drowart, Exsteen and Verhaegen (1964) [3], mass spectrometry. 

^Engelke, Sander and Zare(1976).'[9], chemiluminescent studies. 


PRECEDING page rank not filmed 



- 5 - 


121 Irvin and Dagdiagian(l980) [10], chemiluminescent studies. 

n The value given is the reinterpreted one (see text). Original value was 
D o ^4.67±0.15 eV. 



Table XII. Theoretical spectroscopic parameters (r e , w e ) for the 
A 1 !! state of BeO and the a 3 II states of MgO, CaO, SrO and BaO 


Molecule/ state 


r e (A) a 



w e (cm 

■) a 


SCF 

CI(SD) 

Expt^ 

SCF 

CI(SD) 

Expt^ 

BeO( 1 II) 

1.456 

1.456 

1.463 

1180 

1169 

1144.24 

M g o( 3 n) 

1.852 

1.857 

1- 

00 

• 

T-i 

690 

685 

650 

CaO( 3 n) 

2.105 

2.087 

2.09 9 

561 

542 

556 

SrO( 3 n) 

2.236 

2.217 

2.196 

453 

491 

464 

BaO( 3 II) 

2.347 


2.289 

424 

— 

445.4 


Spectroscopic parameters are based on a parabolic fit to the lowest energy 
points on a 0.1 Bohr grid. 

^Values taken from Huber and Herzberg [69]. 
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Table XIII. Theoretical dissociation energies (D 0 ) and 2 II— 2 E + energy 
separations for the alkali sulfides Li-Rb 


Molecule 

D 0 (eV) a ' b 

T e (cm *) for A 2 S 

LiS 

3.19(3.22) 

5107(5035) 

NaS 

2.50(2.66) 

3814(4033) 

KS 

2.51(2.68) 

1717(1827) 

RbS 

2.39(2.58) 

1240(1280) 


a Dissociation energies are for the X 2 II states with respect to ground state atoms. 
Zero-point corrections based on theoretical u> e values in Table XIV. 

^The SCF values are given first with CI(SD) values given in parentheses. 
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Table XIV. Theoretical spectroscopic parameters (r e , 0J e ) for the 
X 2 n and A 2 S + states of the alkali sulfides Li-Rb a 


Molecule 

r e l 

(A ) b 

“>e( cm 

-> b 


x 2 n 

a 2 e + 

x 2 n 

A 2 S + 

LiS 

2.166(2.163) 

2.100(2.093) 

582(572) 

538(536) 

NaS 

2.515(2.487) 

2.454(2.424) 

397(339) 

326(319) 

KS 

2.890(2.843) 

2.776(2.717) 

223(239) 

250(261) 

RbS 

3.036(2.983) 

2.903(2.834) 

208(214) 

198(200) 


Spectroscopic parameters are based on a parabolic fit to the lowest energy 
points on a 0.1 Bohr grid. 

^The SCF values are given first with the CI(SD) values given in parentheses. 


'i MiAjS 
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Table XV. Theoretical spectroscopic parameters (r ei w e , D e ) for the 
A 1 !! states of BeS, MgS and CaS. a 


Molecule 

r e (A ) b 

w e (cm x ) b 

D e (eV) b,c 

a 3 II-A 1 n energy 
separation(cm _1 ) b 

BeS 

1.926(1.921)^ 

796(784) d 

2.04(2.29) 

937(941) 

MgS 

2.335(2.310) 

431(415) 

1.49(1.70) 

731(859) 

CaS 

2.643(2.601) 

320(338) 

2.24(2.41) 

319(374) 


a The a 3 II r e and u/ e values are not reported since they are essentially the same as 
those for the A 1 !! states. 

^The SCF values are given first with the CI(SD) values in parentheses. 

c Dissociation energies for the A 1 !! states (no zero-point correction) with respect 
to ground state atoms. To convert to ground state D 0 , add TjA^) and sub- 
tract X x E + zero-point energy. 

^Experimental values are r e =1.907 s and u> e =762.46 cm -1 from observation of the 
electronic spectrum by Cheetham et al [77]. 



Table XVI. Dissociation energies (D 0 ) of the alkali 
and alkaline- earth hydroxides 


Molecule 

Theoretical D 0 (eV) 

LiOH 

4.65(4.64) 

NaOH 

3.48(3.51) 

KOH 

3.49(3.59) 

RbOH 

3.43(3.54) 

CsOH 

3.57(3.71) 

BeOH 

4.64(4. 70) e 

MgOH 

3.33(3.31) 

CaOH 

4.10(4.13) 

SrOH 

4.11(4.16) 

BaOE 

4.30(4.38) 


Experimental D 0 (eV) 

4.55±0.09 b , 4.64±0.09 c , 4.54 d 
3.47±0.09 b , 3.38 d 
3.69±0.09 b , 3.53 d 
3.77±0.09 b , 3.61 d 
3.90±0.09 b , 3.95 d 
5.51 f , 4.94±0.43 g 

4.16±0.13 h , 2.43±0.22 1 , 3.59±0.22-j 
3.19±0.22 k 

4.44±0.07 1 , 4.51 m , 4.08±0.13 n 
4.00±0.17° 

4.31±0.10 p , 4.47 m , 4.03±0.13 q 
4.00±0.17° 

4.64 r , 4.94 m , 4.73±0.13 n 
4.55±0.17° 


SCF values are given first with CI(SD) values in parenthesis. We used the fol- 
lowing zero-point corrections (eV) in converting our D e values to D 0 : 
LiOH(O.lO), NaOH(0.08), KOH(0.07), RbOH(0.06), CsOH(0.06), BeOH(O.lO), 
MgOH(0.07), CaOH(0.08), SrOH(0.08), and BaOH(0.08). 

^Cotton and Jenkins(1969) [82], atomic absorption spectroscopy in flames . 

c Zeegers and Alkemade(l970) [83], flame photometry. 

^Jensen(1970) [84], flame photometry. 

e Calculation at linear geometry. Calculations with a smaller basis set show a 
bent geormetry to be about 0.01 eV more stable, 
f 

Ko, Greenbaum and Farber(1967) [85], molecular flow effusion. 

^Inami and Ju(1968) [86], flame photometry, as reported in JANAF 1978. 
^Cotton and Jenkins( 1969a) [88], flame photometry. 
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^ulewicz and Sugden(1959) [89j, flame photometry. 

^Estimate from JANAF Tables(1978) [87]. 

Murad(1980) [90], mass spectrometry. 

^Kalff and Alkemade(1973) [96], flame photometry. 

m Cotton and Jenkins(1968) [92], atomic absorption spectroscopy in flames. 
n Ryabova, Khitrov and Gurvich(l972) [93], flame photometry. 

°Murad(1981) [94], mass spectrometry. 

p Hurk, Hollander, and Alkemade(1974) [95], recalculation of work in footnote 1. 
q Gurvich, Ryabova, Khitrov, and Starovoitov(l971) [91], flame photometry. 
r Stafford and Berkowitz(1964) [97], mass spectrometry. 
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Table XVII. Theoretical spectroscopic parameters (r e , w e ) 
for the ground states of the alkali and alkaline- earth hydroxides 


Molecule 


r e (A) b 



w e( cm 

-) b 


SCF 

CI(SD) 

Expt 

SCF 

CI(SD) 

Expt 

LiOH 

1.576 

1.573 

1.58 c 

941 

931 


NaOH 

1.940 

1.932 

1.95 d 

2.196 f 

588 

579 

431 e 

KOH 

2.235 

2.208 

455 

467 

408 g 

RbOH 

2.349 

2.323 

2.301 n 

380 

398 

354 e 

CsOH 

2.448 

2.419 

2.391 n 

350 

378 

400±80 1 , 336^ 

BeOBi 

1.371 

1.372 j 

— 

1340 

1343 

— 

MgOH 

1.757 

1.758 

1.976 1 

810 

802 

62 3 1 

CaOH 

2.006 

1.988 

612 

629 

SrOH 

2.134 

2.117 

2.10 m 

547 

533 

522 m 

BaOH 

2.251 

2.234 

.... 

483 

502 



a The ground state is 1 £ + for the alkali hydroxides and 2 £ + for the alkaline-earth 
hydroxides. 

^The metal-oxygen values assuming a rigid OH subunit with a fixed OH bond 

o 

distance of 0.94 72 A . The parameters are based on a parabolic fit to the lowest 
energy points on at least a 0.1 Bohr grid. 

c Chase et. al. (1974) [78], computed using moment of inertia reported by Freund 
et. al., 25th Spectroscopy conference (Columbus, 1970) paper E8. 

^Kuijpers, Torring and Dymanus(1976) [101], microwave. 

e Acquista and Abramowitz(1969) [102], matrix isolation. 

f 

Pearson, Winnewisser and Trueblood(1976) [103]. 

g Belyaeva, Dvorkin and Sheherba(1966) [104]. 

^Lide and Matsumura(1969) [105]. 

^ide and Kuczkowski(1967) [106], gas phase. 

■^Acquista, Abramowitz and Lide(1968) [107], matrix isolation, 
k 

Computed at the linear geometry. Calculations with a smaller basis yield a 
bent structure at the SCF(CI) level of R(Be-O)=1.390A (1.396A ), R(0- 
H) =0.933 A (0.948 A ) , and 0=34 ‘ (33 ° ) . 

^Hilborn, Qingshi, and Harris(1983) [108], gas phase. 
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m "Nakagawa, 


Wortnsbecher 





Table XVIII. Theoretical dissociation energies for some 
alkaline- earth (iso)cyanides 


Molecule 


dissociation energies (eV) a 



Cyanide 

Bmp° 

Isocyanide 

Expt c 

BeCN 

3.62(3.83) 

3.17(3.38) 

4.10(4.09) 

4.46±0.29 

MgCN 

3.07(3.22) 

3.00(3.12) 

3.36(3.35) 

4.21±0.26 

CaCN 

3.78(3.89) 

3.88(3.96) 

4.12(4.09) 

4.46±0.22 

BaCN 

4.02( ) 

4.17( ) 

4.34(4.33) 

4.87±0.26 


a Dissociation energies (without zero-point corrections) relative to CN( 2 E + ) + 
M( 1 S). SCF values are given first with CI(SD) values in parentheses. 

^The metal atom is approaching the CN bond mid point. 

c L’vov and Pelieva(1980) [118], electrothermal atomic absorption spectrometry. 


v 
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Table XIX. Recommended dissociation energies (D 0 ) for selected alkali 
and alkaline-earth fluorides, chlorides, oxides, sulfides, hydroxides 

and isocyanides. 


Molecule 

State 

D 0 (eV) 

Molecule 

State 

D 0 (eV) 

LiF 

*E + 

6.06 

BeF 

2 S + 

5.94 

NaF 

l E + 

5.00 

MgF 

2 E + 

4.68 

KF 

*E + 

5.10 

CaF 

2 S + 

5.53 

RbF 


5.07 

SrF 

2 E + 

5.62 

CsF 

1 E + 

5.27 




LiCl 

X E + 

4.89 

Bed 

2 e + 

3.90 

NaCl 

1 E + 

4.22 

Mgd 

2 e + 

3.26 

KC1 

*E + 

4.37 

Cad 

2 e + 

4.14 

RbCl 

*E + 

4.36 

Srd 

2 e + 

4.23 

CsCl 

*e + 

4.58 




LiO 

2 n 

3.84 

BeO 

a e+ 

4.69 

NaO 

2 n 

2.83 

MgO 

1 E + 

2.75 

KO 

2 e + 

2.86 

CaO 

*E + 

4.14 

RbO 

2 S + 

2.90 

SrO 

i s + 

4.36 

CsO 

2 e + 

3.10 

BaO 

*E + 

5.75 

LiS 

2 n 

3.30 

BeS a 


2.29 

NaS 

2 n 

2.67 

MgS a 


1.70 

KS 

2 n 

2.74 

CaS a 

l n 

2.41 

RbS 

2 n 

2.66 




LiOH 

l E + 

4.64 

BeOH 

2 e + 

4.70 

NaOH 

l E + 

3.51 

MgOH 

2 S+ 

3.31 

KOH 

l E + 

3.61 

CaOH 

2 S+ 

4.15 

RbOH 

*e + 

3.58 

SrOH 

2 e+ 

4.18 

CsOH 

*E + 

3.77 

BaOH b 

2 S+ 

4.60 




BeNC 

2 E+ 

4.15 




MgNC 

2 E+ 

3.37 
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CaNC 2 E + 4.13 

BaNC b 2 E + 4.50 


The dissociation energy (D e ) reported is for the excited A IT state and does not 
include a zero-point correction (see Table XV). 

^The results are less accurate (~0.2 eV uncertainty) because the core electrons 
are described by a relativistic effective core potential (see text). 



1 - 


Figure caption 

Figure 1. The dissociation energy (without zero point corrections) with respect to the 
ionic limits versus bondlength. 
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